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Abstract 

We continue our investigation of versality for parametrized families of linear retarded 
^ ' functional differential equations (RFDEs) projected onto finite-dimensional invariant man- 

ifolds. In this paper, we consider RFDEs equivariant with respect to the action of a com- 
pact Lie group. In a previous paper (Buono and LeBlanc, to appear in J. Diff. Eqs.), 
we have studied this question in the general case (i.e. no a priori restrictions on the 
RFDE). When studying the question of versality in the equivariant context, it is natural 
to want to restrict the range of possible unfoldings to include only those which share the 
same symmetries as the original RFDE, and so our previous results do not immediately 
apply. In this paper, we show that with appropriate projections, our previous results 
on versal unfoldings of linear RFDEs can be adapted to the case of linear equivariant 
RFDEs. We illustrate our theory by studying the linear equivariant unfoldings at double 
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Hopf bifurcation points in a Ds-equivariant network of coupled identical neurons modeled 
by delay-differential equations due to delays in the internal dynamics and coupling. 



1 Introduction 

Symmetry plays an important role in the description of many physical systems. Whether 
the symmetry is an intrinsic part of the physical system, or whether it is merely a modeling 
assumption, the mathematical model for any such system must reflect this symmetry. When 
the mathematical model is a differential equation, this is achieved by imposing equivariance 
of the differential equation. That is, if F is a group of transformations which represents the 
symmetries of the physical system, then the differential equation model must commute with an 
action of F. An important consequence of F-equivariance is that any solution of the differential 
equation is mapped onto other solutions by the elements of F. Thus, the algebraic structure of 
the group F yields important information about the geometrical structure of solutions of the 
differential equation in phase space. 

The theory of equivariant dynamical systems is by now well-developed and understood (see, 
for example [6] and [7]). F-equivariance imposes strong algebraic restrictions on the functional 
form of the differential equation. This has important implications, for example, in the study 
of local stability and bifurcation of equilibrium points. In particular, since the linearization of 
the vector field at the equilibrium point must commute with F, then this linearization usually 
has eigenvalues of high geometric multiplicity. 

In this paper, we will be interested in systems which are modeled by parametrized families 
of retarded functional differential equations 



where C{a) is a parametrized family of bounded linear functional operators and {u, a) i— > 
F{u, a) is nonlinear and smooth with F{0, a) = and D„F(0, a) = for all a, such that (1.1) 
is equivariant under the action of a compact Lie group of symmetries F for all a. Examples of 
such systems can be found in the study of neural networks modeled by symmetrically coupled 
cell systems with delays in the couphng ([11, 12, 14, 15]), in the study of semiconductor lasers 
subject to optical feedback [10, 13] or symmetrically coupled arrays of lasers [9]. Specifically, 
we arc interested in the following question regarding (1.1): consider the equilibrium a; = 0, 
a = of (1.1), and let y be a finite-dimensional invariant manifold for the dynamics of (1.1) 
which passes through this equilibrium point. Let 



z{t) = C{a) zt + F{zt, a), ae , 



(1.1) 



X — B{a)x -\- G{x, a) 



(1.2) 
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be the parametrized family of ordinary differential equations whicfi gives the dynamics of (1.1) 
restricted to V, where B{a) is a parametrized family of matrices, and G is the nonlinearity. 
Under what conditions on C{a) in (1.1) can we be assured that B{a) is a versal unfolding of 
the matrix B{0) within the space of T-equivariant matrices! Recall that, loosely speaking, the 
family B{a) is a versal unfolding of i?(0) if it is, up to similarity transformations and changes of 
parameters, the most general C°° perturbation of B{0) (a precise definition is given in [2] in the 
absence of F-equivariance, and in Section 3 of this paper in the F-equivariant case). The above 
question is especially important in the case where (1.2) represents a center manifold reduction 
of (1.1). If B{a) is not versal, then this could lead to restrictions on the bifurcation behavior 
of (1.1) near the equilibrium. It is clear that answering this question does not depend on the 
nonlinearity F in (1.1); therefore, we restrict our attention to the class of parametrized families 
of F-equivariant linear retarded functional differential equations. 

This problem was studied in [2] in the absence of F-equivariance. In particular. Theorem 
7.4 of that paper gives a method for constructing a parametrized family C{a) of bounded 
linear functional operators such that the family of matrices B{a.) in (1.2) is a versal unfolding 
of the matrix B{0). Our approach in this paper will be to construct a suitable equivariant 
projection of the family C[a) of Theorem 7.4 of [2], and to show that the resulting associated 
F-equivariant family of matrices B[a) in (1.2) is a versal unfolding of B[0) within the space of 
all F-equivariant matrices. 

The paper is organized as follows. In Section 2 we review some basic theory of linear 
retarded functional equations while Section 3 contains the results concerning versal and mini- 
versal unfoldings of equivariant matrices. In Section 4, our main result about the construction of 
equivariant unfoldings for linear retarded functional differential equations is presented. Finally, 
in Section 5 a detailed computation of the equivariant mini- versal unfolding of a Da-equivariant 
linear delay-differential equation at a double Hopf bifurcation is discussed in both cases of 
simple and double imaginary eigenvalues. 

It is assumed that the reader is familiar with the terminology and results of [2] , [3] , [4] , [5] 
and [8]. 

2 Reduction of equivariant linear RFDEs 

Let Cn = C{[—T, 0], C") be the Banach space of continuous functions from the interval [— r, 0], 
into (r > 0) endowed with the norm of uniform convergence. We will first consider the 
linear homogeneous RFDE 

z{t)^Co{zt), (2.1) 
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where Cq is a bounded hnear operator from C„ into C". We write 

Co{^) = f drj{9Me), 

where 77 is an n x n matrix- valued function of bounded variation defined on [— r, 0] . 

We will suppose that (2.1) has certain symmetry properties. Specifically, we let F be a 
compact group of transformations acting linearly on C". We say that (2.1) is r-equivariantii 

7 ■ dr]{e) = dr]{e) -7, V7 G r, ^ G [-r, 0]. (2.2) 

A consequence of this definition is that if x{t) is a solution to (2.1), then so is 7 • x{t), for all 

7 G r. 

Suppose that A C C is a non-empty finite set of eigenvalues of the infinitesimal generator 
Aq for the semi-fiow of (2.1), i.e. A G A if and only if A satisfies the characteristic equation 

/■O 

detA(A)=0, where A(A) =A4-y dr]{e)e^\ (2.3) 

where /„ is the n x n identity matrix. Using adjoint theory, it is known that we can write 

Cn^PeQ (2.4) 

where the generalized eigenspace P corresponding to A and the complementary subspace Q are 
both invariant under the semifiow of (2.1), and invariant under Aq. 

Define C* = C([0, r], C"*), where C^* is the n-dimensional space of row vectors. We have 
the adjoint bilinear form on C* x C„: 



(^, = mm -J I ^(^ - o)d7^{eMOd^. (2.5) 

We let $ — ((/?!, . . . , (fc) be a basis for P, and ^ — co^V^i, . . . jipc) be a basis for the dual 
space P* in C*, chosen so that (\1/, $)„ is the c x c identity matrix, !(.. In this case, we have 
Q = {ip & Cn '■ i"^, ^p)n = 0}. It follows that $' = $5 where i? is a c x c constant matrix. The 
spectrum of B coincides with A. Using the decomposition (2.4), any 2; G C„ can be written as 
z = ^ X + y, where x G C'^ and y E Q is a. function. The dynamics of (2.1) on P are then 
given by the c-dimensional linear ordinary differential equation 

X = Bx. (2.6) 

The important point, as we will now show, is that F-cquivariancc of (2.1) forces F-invariance 
of the splitting (2.4) and F-equivariance of the reduced ordinary differential equation (2.6). 
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Proposition 2.1 The spaces P, Q and P* defined above are V invariant, and the cx c matrix 
B in (2.6) commutes with a representation G : F — > GL{c, C) of the group F. 

Proof We start by noting the following identity which follows trivially from F-equivariance 
of (2.1), and from (2.5): 

(*, 7 •$)„ = (*• 7, $)„ V7eF. (2.7) 

The functions ipi,...,(fc which span P are generalized eigenf unctions of the operator A 
defined by Aip — (f' whose domain is 

V{A) = I 99 e C„ : e C„ and ^^'(0) = J dr){e) if {6) | . 

Let 7 G F, (/? G P, and consider = 7 ■ Then obviously (p' G C„ since ip' G C„. Moreover, 
from (2.2) and the fact that G I^(^), we get 

^'(0) = 7 • ¥^'(0) = 7 r ^^(^) ^(^) = f dr0) 7 ■ m = f dv{9) m, 



— T 



SO (fi E ^{A). Now, there exists an integer k and a A G A such that {A — \I)^ip ~ 0. Since 
A{'^-ip) — (7 • (^)' = 7 • = 7 • A (/?, it is easy to see that (A — XI)^(p — 0, and we conclude that 
(f E P. Thus, P is F-invariant. A similar argument shows that P* is also F-invariant. This 
implies that for all 7 G F, there exist c x c matrices ^(7) and H{'y) such that 

7 . $ = $ . G(7) and * • 7 = H{^) ■ (2.8) 
In fact, it follows from (2.7) and the fact that (^', $)„ = Ic that 

i7(7) = i?(7)(*, = ■ $)„ = (* • 7, = (*, 7 • = 

(vl/,$.G'(7))n = (^,$)nG'(7) = G'(7). 

Note that 

Giaj) = (*, $)„G'((77) = (*,$• G'(a7))„ = (*, (77 • $)„ = 
(^ • a, 7 • *)n = (G(ct) • ^, $ • G(7))n = G(^)(^, *)nG(7) = G{a)G{^), 
so that the mapping 

F — > GL{cX) 
7 ^ G(7) 
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is a representation of the group T in the space of invertible cxc matrices. 
Now, (fi e Q ii and only if — 0. Then using (2.7), one gets that 

7 ■ = ■ 7, = (^(7) ■ ^, ip)n = G'(7)(^, = 0, 

so Q is F-invariant. 

FinaUy, it foUows from the fact that B = {'^, $')„ that 

B ■ G{^) ^ G{^) ■ B V7er. (2.9) 



3 Equivariant versal unfoldings of matrices 

Let B e Matcxc, where Matcxc denotes the space of c x c matrices with entries in C. Recall 
that a p-parameter unfolding of S is a C°° mapping B : — > Matcxc such that B(q;o) = B 
for some e C^. We then say that a p-parameter unfolding B{a) of i? is a versal unfolding of 
B if, given any ^-parameter unfolding A{P) of B (with A{Po) — B), there exists a C°° mapping 

: — >CP 

with 0(/9o) = cto, and a g-parameter unfolding C{(3) of the identity matrix such that 

A{(3) = C{P)B{4>{P)){C{P))-\ V/3 near /3o. 

If ;B is a versal unfolding which depends on the least number of parameters, then B is called 
a mini-versal unfolding. Thus, one may view a versal unfolding of B as the most general C°° 
perturbation of B up to similarity and change of parameters. One then has the following 
sufficient criterion [1] for establishing the versality of a given unfolding of B: 

Proposition 3.1 Let B{a) be a p-parameter unfolding of B & Matcxc- Let E denote the 
similarity orbit of B, i.e. 

E = { CBC-^ : C e GL{c, C) } C Matcxc- 

// 

Matcxc = TbE + D^B{ao) ■ (3.1) 

where TbT: is the tangent space to E at B, then B is a versal unfolding of B. If in addition, 
the codimension ofTs^ in Matcxc equal to p, then B is mini-versal. 
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In this paper, we are not interested in the most general perturbation of B, since we are 
only considering matrices in Mat^xc which commute with the group F. In fact, we need only 
consider perturbations in the subspace Mat^^^ defined by 

Mat^^, = { A e Matexc : ^(7)^ = ^^(7), V7 e T }, 

where G is the representation of F in GL{c, C) defined in (2.8). We therefore adopt the following 
definition: 

Definition 3.2 A p-parameter unfolding B{a) of B G Mat^^^ is said to be a F-unfolding of B 
if B{a) G Mat^^^ for all a. The T-unfolding is said to be F-vcrsal if, given any q-parameter 
F-unfolding A{P) of B (with A{(3o) = B), there exists a C°° mapping 

: C — ^ 

with (f){Po) = cto, (ind a q-parameter F-unfolding C(/3) of the identity matrix such that 

A{P) = C(/3)-B(0(/3))(C(/3))-\ V/? near Po- 

If B is a F-versal unfolding which depends on the least number of parameters, then B is called 
a F-mini-versal unfolding. 

The following lemma, which is the F-equivariant version of Proposition 3.1, is proved in Ap- 
pendix A. 

Lemma 3.3 Let B{a) be a p-parameter F-unfolding of B G Matexc- denote the F- 

similarity orbit of B, i.e. 

TF ^{ CBC-^ : C e GL{c, C) n Mat^^, } C Mat^^,. 

// 

Mat^e = TbTF + D^B{ao) ■ (3.2) 

where TbTF is the tangent space to at B, then B is a F-versal unfolding of B. If, in addition, 
the codimension ofTs^F in Mat^^^ is equal to p, then B is F-mini-versal. 

Lemma 3.4 There exists a projection : Matexc — ^ ^^'^^xc ■^^^^ ^^^^ 

7rl(AM)^A7rl(M) and nl(MA) ^ Trl(M) A, VA e Mat^^,, M e Matexc- (3.3) 
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Proof For M e Matcxc, define 



where the integral denotes the Haar integral on F, which we assume has been normalized so 
that the measure of F is 1 (recall that we assume that F is compact). Since G is a morphism 
and since the Haar integral is translation invariant, we have that for all u e F 

G{a)i:l{M){G{a))-' = ^ G{a^)MG((a^)-') = 7rl(M), 

so that 7rf (M) G MatJ^^c- Moreover, we clearly have that 7rf (M) = M if and only if M G Matly.^. 
We then conclude that Trf is idempotent, with range Mat^^^. Finally, (3.3) follows from the 
definition of tt^ and the fact that A commutes with G(7), for all 7 e F. I 



Proposition 3.5 Let B be a p-parameter unfolding of B E Mat^^c such that (3.1) holds (which 
implies by Proposition 3.1 that B is a versal unfolding of B). Let tt^ be the projection 0/ Matcxc 
onto Matcxc such as in Lemma 3.4. Then B^ = tt^ (S) is a F-versal unfolding of B. 

Proof We first show that T^ciT^s) = TTF^. Let x = [B,y] e TS^, where y G Mat^^c- 
Since obviously x G Matcxc, ^^^^ ^ — '^ci^) — '^cH^tV]) ^ tTcC^^b), since [B,y] is also an 
element of TE^. Thus TS^ C 7t^{TEb). Now let x G 7rJ'(rEB), i.e. x = 7r^{By-yB) for some 
y G Matcxc- Then, from Lemma 3.4, we have 

X = TrliBy - yB) = i:l{By) - 7^l{yB) = [B^M] ^ TE^. 



Therefore, we have shown that 7r^(TEs) 

Now, we will show that 7r^(L>«'B(ao) -C^) = D^B^{ao)-CP. We have that x G DaB^{ao) ■ C^, 
if and only if there exists v E such that 



X 



^^{B^iao + ev)) 



£=0 



-^{7rl{B{ao + ev)) 



e=Q 



de 



{B{ao + ev)) 



£=0 



G7rr(L'«B(ao)-n. 



Therefore, 7r^{D„B{ao) ■ CP) = DaB^{ao) ■ CP. 
Finally, we show that 



Mat! 



TEl + DaB'^{ao)-CP. 
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Obviously, from the preceding results, we have that 

Let X G Mat^x^, then from (3.1) we have that x — Xa + xj,, where Xa G TE^ and xi, G 
D^B{ao) ■ CP. Therefore 

X = nlix) = TT^ixa) + TT^ix,) G ttHTEb) + T:l{D^B{a^) ■ CP) = TE^ + D^B^ia^) ■ CP. 

Thus, we have shown that Mat^^c — + D^B^ {ao) ■ C^, and by Lemma 3.3, we conclude 
that B^ is a F-versal unfolding of S. ■ 

Remcirk 3.6 Let B G Mat^^o given, and let B{a) he a p-parameter unfolding of B which 
is mini-versal. Then (3.1) holds as a direct sum, and by the previous Proposition, B^{a) = 
nl{B{a)) is a T-versal unfolding of B. It should be noted that in general, B^{a) need not be 
r -mini-versal, even though B is mini-versal (in fact it is quite easy to constrct examples to 
illustrate this fact). However, it is always possible to extract a V-mini-versal unfolding of B 
from B^ , using a straighforward procedure as we now illustrate. 

Let W be any direct sum complement of TbS in Matcxc, i-e. 

Matexc = Tb^ © W. (3.4) 

For example, if W is chosen to be the space which is spanned by the matrices Qi,...,ftp 
described in section 6 of [2], then the decomposition of any cxc matrix using (3.4) is particularly 
simple to compute. Now, for each z = 1, we define 

S[^^(i3>o + ee.))U, 

where Cj is the p-dimensional vector whose i*'* component is 1 and all other components are 
zero. Note that 

D^B^{ao)-CP^C-{Bl,...,Bl}. 
Now, decompose these matrices using (3.4), i.e. write 

p 

= i^l,...,p, (3.5) 

where yi G Matcxcj and the 9ij are scalars. Let © = {9ij) be the p x p matrix coming from 
the decomposition (3.5). If we then let rows ii, . . .ik of the matrix 6 determine a maximal set 
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of linearly independent vectors in the rowspace of the matrix ©, it follows from simple linear 
algebra that 

MatL, = TBE^©C-{4';,...,5[J, 
and so the following /c-parameter F-unfolding of B is F-mini-versal: 

Recall that the RFDE (2.1) is defined on C", and that F is a subgroup of GL{n, C). It will 
be useful to define a projection analogous to Trf on Mat„xn- 

Definition 3.7 Let A e Mat^xn- We define 

Obviously, it follows that for all A e Mat^xn; we have a7r^{A) — n^{A) a for all cr e F. 
Proposition 3.8 Suppose C e Matcxc is such that 

C = *(0)A$(t*), (3.6) 
where A is some n x n matrix, and r* is some fixed number in [— r, 0] . Then 

i:l{C) = ^{Q)i^'^{A)^{T*). 

Proof This is just a simple computation, using (2.8): 

7r[(C) = j G'(7)^(0)A$(r*)G(7-^)d7 = ^(0) {^j 7/17-^^7^ $(r*). 
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4 Construction of a F-equi variant A-versal unfolding 

Consider now a C°° parametrized family of linear RFDEs of the form 

(4-1) 

where a e C^, and C{ao) = Cq is as in (2.1), for some ag e C^. In the sequel, we will assume 
that a translation has been performed in the parameter space such that ao — 0. It was 
shown in [2] that for a given set A of solutions to (2.3), a parameter-dependent reduction of 
(4.1) to the c-dimensional generalized eigenspace corresponding to A is given by 

X = B(a)x, (4.2) 

where 

B{a) = B + ^(0)[£(a) - £o]($ + h{a)), (4.3) 
where h is a smooth n x c matrix- valued function such that h{0) = 0, and B is as in (2.6). 

Remark 4.1 The function h is the unique solution (via the implicit function theorem) of a 
certain nonlinear equation N{a,h) = (see equation (4-5) in Proposition 4..! of [2]). It is a 
straightforward computation that for all h and all ^ eT, we have 

N{a, 7 h G(7-^)) = 7 N{a, h) G{^-^) 

so that if h{a) is such that N{a, h{a)) — 0, then by uniqueness of the solution, we must have 
^ h{a) G(7~^) = h{a), i.e. 

for all a near and for all ^ eT. Thus, if (4-1) is V-equivariant for all a, then it follows from 
(4-3) that B{a) e Mat^^^^ for all a near in C^, i.e. B is a T-unfolding of B. 

We recall the following definition from [2]. 

Definition 4.2 The parametrized family of RFDEs (4-1) is said to be a A-versal unfolding 

(respectively a A-mini-vcrsal unfolding^ for the RFDE (2.1) if the matrix B{a) defined by (4-3) 
is a versal unfolding (respectively a mini-versal unfolding) for B. 

In the case where (2.1) is F-equivariant, we have 

Definition 4.3 The parametrized family (4-i) is said to be a F-equivariant A-versal unfolding 
(respectively a F-equivariant A-mini- versal unfolding^ for (2.1) if (4-1) is F-equivariant for all 
a, and the matrix B{a) defined by (4-3) is a V-versal unfolding (respectively a V-mini-versal 
unfolding) for B. 
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In [2], we showed that given B e Matcxc, there exists points ro,ri, . . . ,Tc^q in [— r, 0], and 
nxn matrices A^, j = 0, . . . ,c — q, m = 0, . . . ,S (where S is the codimension of TE^ in Matcxc) 
such that the following 5-parameter unfolding of B 

5 /c-q \ 

B{a) = 5 + E E *(0)A7*(r,) (4.4) 

m=l \i=0 / 

satisfies (3.1) as a direct sum (i.e. i3 is a mini-versal unfolding of B). From this, it follows (see 
Theorem 7.4 of [2]) that if C{a) is the 5-parameter family of bounded linear operators from 
C([— r, 0], C") into C" defined in the following way: 

S /c-q \ 

m=l \j=0 J 

where the are complex parameters and Cq is as in (2.1), then (4.1) is a A-versal unfolding 
of (2.1). 

We now state and prove the main result of this paper. 

Theorem 4.4 Suppose B G Mat^^^, and consider the matrices AJ^ which appear in the versal 
unfolding (4-4) of B. Then the family jC{a) of T-equivariant bounded linear operators defined 
by 

5 /c-q \ 

£{a) z = £oz+Y,^m(Yl <i^T) 'i^j) ' (4-5) 

m=l \i=0 / 

(where vr^ is as in Definition 3.7) is such that (4-1) is a T-equivariant A-versal unfolding of 
(2.1). 

Proof From Propositions 3.5 and 3.8, we conclude that the versal unfolding (4.4) of B is 
such that ^ 

m=l \j=0 

is a F- versal unfolding of B. Thus, if we define C{(y) as in (4.5), then B[a) (as in (4.3)) is such 
that 

B{a) = Til{B{a)) + ^(0)[£(a) - Co\{h{a)). 

We conclude that B{a) is also a F-versal unfolding of 5, since /i(0) = implies that DaB{f)) = 
D^7il{B{a))\a=o. ■ 
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5 An example with F = B3 



In this section we compute the versal unfolding of a Da-equivariant system of hnear delay- 
differential equations near points of double Hopf bifurcations. Consider the group D3 generated 
by {k, 7} with — 'j^ — 1 and the representation p : D3 — > GL(3, C) given by 



p{k) =001, p(7) =001. (5.1) 





Wu et al. [15] and Ncube et al. [12] consider the Da-equivariant linear delay-differential 
equation 

= -Uj{t)^aUj{t-Ts)+ (5[Uj-^i{t-Tn)+Uj+i{t-Tn)\, j ^ 1,2,?, (5.2) 

where the indices are taken mod 3. This linear system appears as the linearization at an 
equilibrium of a coupled system of three multiple-delayed identical neurons. In those papers, 
the focus is on finding periodic solutions via the Hopf bifurcation theorem. 

The spectrum of the linear operator associated to (5.2) is obtained by solving the charac- 
teristic equation 

5(A) = Ai(A)A2(A)2 = 0, 

where 

Ai(A) = A + 1 - ae-^^^ - 2^6"^^", 
A2(A) = A+ 1 -ae-^^» +/?e-^^". 

We now sketch an argument showing that points of double Hopf bifurcation exist for both 
factors Al and A2 of the characteristic equation. 

Suppose that we are looking for a point of double Hopf bifurcation solution to Ai = 0, then 
the eigenvalues iu) must satisfy 

a; -|- Q;sin(a;rs) = — 2/?sin(a;r„) 
1 — Q;cos(a;T5) = 2(3 cos{ujTn) . 

It is easy to solve these equations for a and in terms of the other parameters, we obtain the 
following two equations 

a = ±^(1 - 2(3 cos{ujTn)Y + (w + 2/3sin(a;r„))^ 

1 f —uj — 2(3sniiujTn) 

~ — arctan ' 



uj V 1 ~ 2(3cos{ujTn) 
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Figure 1: Curves of imaginary eigenvalues in (a, r^) -space for (a) Ai = where (3 — —0.5 and 
Tn — 4, and (b) A2 = where P = 0.5 and r„ = 3. 

Then a and with (3 and r„ fixed represent curves parametrized by co in {a, rs)-space. Note that 
there are several branches of curves depending on the branch of the arctan which is chosen. 
Figure 1(a) shows some of these curves for P — —0.5 and t„ = 4. Points of double Hopf 
bifurcation lie at the intersection of these curves. Points of double Hopf bifurcation for A2 = 
are found in a similar way by isolating a and Ts as functions of /3, Tn and u. Figure 1(b) shows 
several curves of imaginary eigenvalues plotted for P = 0.5 and r„ = 3. 

From the graphics in Figure 1, it is reasonable to suppose that the characteristic equation 
has a point of double Hopf bifurcation at {a* , P* , t* , r*) with r*, r* nonzero, r* 7^ r* and 
frequencies ui 7^ uj2- Of course, the point {a* , P* , r* , r*) as well as the frequencies are different 
whether Ai = or A2 = 0, however, to keep the notation to a minimum we use this unique 
label for both cases. 

We let A = {±iLOi, ±iuj2} with lvi ^ and a;2 7^ be the set of eigenvalues on the imaginary 
axis at a double Hopf bifurcation point. If the double Hopf point occurs for the Ai factor, 
then the eigenvalues of A are simple and the restriction of the linear operator to the center 
eigenspace is given by 

/ Uli \ 
-cuii 
~ a;2i 
\ -cu2i J 

If the double Hopf point occurs for the A2 factor, then the eigenvalues are double and the hnear 
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flow on the center eigenspace is given by the following matrix 



/ 


ujii 




















\ 







cuii 




























—cuii 




























—Uli 

























































































\ 























-0021 ) 



We now compute the Ds-equivariant A-mini-versal unfolding of the linear bounded operator in 
both cases. Recall that our goal is to flnd matrices AJ- from which we can write (4.5). The 
computations rely on the results in [2], mainly Section 6 and Lemma 7.1. 



5.1 Simple imaginary eigenvalues. 

The center eigenspace P has basis 



= (we"^i^ Me-"*'!^ ue"^^\ ue-'"^^^) 



where 



u 




Its dual, P*, has basis ^(s) = (V'll'S), '02('S), "03(5), '04(s))* where -02 = '^i, '04 = i>z and i'j{s) 
(V'ji(s), -0^2(5), -0^-3(5)), j = 1,2,3,4. 



Basis of the unfolding space. Let W be the unfolding space spanned by matrices Qi, . . . , 
defined in Section 6 of [2] . The dimension of W for matrix B is given by 

X^f^(2£-l)n,-, (5.3) 

i=i 1=1 

where r = 4 is the number of distinct eigenvalues, kj = 1 is the number of Jordan blocks for 
each eigenvalue and nj,i = 1 is the size of each Jordan block, thus the unfolding space is 4- 
dimensional. The space ^ (W) lying in the range of ^'(O) is linearly isomorphic to W and is also 
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an unfolding space, see Section 6 of [2]. A basis of S{W) is given by '^{0)RP for j = 1,2,3,4 
where R^, R^, R^, i?"^ are 3x4 matrices defined following the procedure established in [2], 
Section 6.3. In [2] we introduce projections Ilj : ^ C {j — 1,2,3,4) with 

Let Vj = e {j = 1,2,3,4) be such that Ilj{vj) = 1 and set R^ = (vi, 0,0,0), 

= (^0, V2, 0, 0), = (0, 0, V3, 0), and i?^ = (0, 0, 0, v^). Note that since -02 = '^i and ^4 = 

we choose V2 — W and V4 — V3. 

Since rank $(0) = 1, by Theorem 7.2 of [2] we need 4 distinct delays tq, ti, T2, T3 to find 3x3 

matrices A*- solving 

3 

i?^ = J]A}$(r,). (5.4) 

j=0 

Projection of the basis. By Theorem 4.4, the solutions to (5.4) are used to write the 
equivariant A-unfolding by projecting the matrices A™ with tt^. Multiplying (5.4) by ^'(O) we 
obtain equation (3.6) with C = \l'(0)i?. Then, projecting C with n^, we obtain 

7tI(C)^ j G(7)*(0)i?G(7~')c^7 = / *(0).p(7)i?G'(7~')c^7 = *(0) J p{i)RG{r')d^. 

(5.5) 

We project the matrices R^ on the subspace of Da-equi variant 3x4 matrices using the projection 
operator obtained from (5.5); that is, since the G representation is trivial on P and F = D3 is 
finite, then the projected matrices are 

Note that the factor | is not necessary to the projection of the matrices and could be removed 
without further consequences. Thus, we obtain new matrices 

















( 


^2 









Vi 








I) 







??2 





:) 


































^ 





















m 







(: 































^4 / 
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where ijj — ^{vj + v'j + v|) and of course r]2 — fji and 774 = 773. Because the matrices R' are 
Ds-equivariant matrices, we may take matrices Aj^ Ds-equivariant with respect to the action of 
D3 on given by (5.1); that is, 

= I of 6f I . (5.6) 





K 


K 








bf 

























Ds-mini-versal unfolding Before solving equation (5.4), we verify which of the four matrices 
R' [j = 1, 2, 3, 4) are necessary to obtain a mini-versal unfolding. We follow the steps described 
after Remark 3.6. Since S{yV) is an unfolding space for B, we set B^' = ^(0)B^ for j = 1, 2, 3, 4. 
The space W is spanned by 4 x 4 matrices (i = 1, 2, 3, 4) where Qi has a 1 on the diagonal 
at {i, i) and is zero elsewhere. Then, it is easy to compute that 

Bf-' = [S, 2/1]+ 771(^11(0) + ^12(0) +Vi3(0))l^i 



5f = [5,7/2] +r72(V'ii(0) + V'i2(0)+Vi3(0))Jl2 

= [S, 7/3] +??3(^3l(0) + 7/^32(0) +^33(0))fi3 

B^' = [B,y4 + rj^{i^ + i^ + ^^f^)n^ 

for some matrices 2/1, 2/2, ?/3, ?/4 G Mat4x4. Since the matrix © (defined using (3.5)) is diagonal, 
the four rows are needed to span the rowspace of ©. Hence, 

Mat^^4 = TbE^^ c . |5»3^ B^^,Bf',Bf'} 

and we need to consider all matrices i?' (j = 1, 2, 3, 4) in order to compute the D3-equi variant 
A-mini-versal unfolding. 

Computation of the matrices AJ^. We now solve for the matrices A^^ using equation (5.4). 
Since 

Ap^T,) = (a™ + 26™)$(r,), 
it is sufficient to solve the first line of (5.4) to obtain the matrices A"^. Consider the first line 



of equation (5.4) for R 



(771,0,0,0) = ^(a], 6), 6])$(r,). 

3=0 



17 



This equation can be rewritten as 






V J 



I 



M 



2hl \ 
2b\ 



(5.7) 



where 



M 



( 



-lOJlTO 



-10)2 TO 



-IWITI 



3ia)2n 



-IWIT2 



-IU)2T2 



^IU!2T3 



-IUI2T3 



(5.8) 



Obviously, det M = if any two of the delays are equal or if ui = U2- Now, since we take our 
four delays to be distinct and it is assumed that cui u}2, then generically M is nonsingular. 



Thus, we obtain Li — Ylj=o'^]^{'''j) by solving (5.7). Now, the equation for yields. 



/_0 \ 

Vi 


V J 



= M 



f al + 2bl\ 

a\ + 2h\ 
al + 2hl 



which if we interchange rows 1 with 2 and 3 with 4 corresponds to the complex conjugate of 
system (5.7). So, we can choose and 6| so that L2 — Li. In a similar manner we obtain L3 
and L4 — L2. 

Looking at equation (5.7), we see that for each j = 0, 1, 2, 3, either or 6j can be chosen 
arbitrarily. We exploit this freedom to choose constants and bj^ which preserve the structure 
of equation (5.2). 



Ds-equivariant A-unfolding. Take Tq = 0, n = r;, T2 = t*, 6™ = 0, 6^ = and = for 
m — 1,2,3,4 (recall that we have assumed t*,t* are nonzero and distinct so that matrix M 
in (5.8) is generically nonsingular with this substitution). Hence, 



Liz ^ allz{0) + a\Iz{T*) + 



L3Z = alIz{Q) + alIz{T;) + 






bl 


bl 


bl 





bl 


bl 


bl 








bl 


bl 


bl 





bl 


bl 


bl 






z{T:) + alIz{Ts) + 



bl bl 



bl 

bl 



bl 

bl 



bl bl 



z{T*) + alIz{Ts)+ \ bl bl 



bl bl 
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and the complex Da-equivariant A-versal unfolding is defined by 

4 

C{a) Z = CqZ+Y^ OimLm 
m=l 

4 4 4 



m=l 



m=l 



b^ 

b^ b^ 
b"^ bV^ 



6^ fe^T 



+ J] a^a™/^(r3) + J] I 6^ 6^ h(^3)- 



m=l 



m=l 



6^ 6^ 



By the form of equation (5.7), we can simplify the unfolding even more by setting for instance 
= for m = 1, 2, 3, 4. Thus a Da-equivariant A-versal unfolding of (5.2) near (a*, /?*, r*, r*) 
is given by 

ui{t) = (-1 + ei)ui{t) + (a* + e2)«i(i - <) + (/5* + e3)(w3(i - r*) + W2(t - r*)) 

+ e4(M3(^ - T3) + W2(t - T3)) 



tt2(i) = (-1 + ei)«2(t) + («* + e2)u2{t - r*) + {13* + e3)(Mi(t - r*) + U:,{t - r*)) 



(5.9) 



W3(t) = (-1 + ei)M3(t) + (a* + e2)u3{t - t*) + (/3* + e3)(w2(t - r*) + ui{t - r*)) 
+ e4(M2(t - T3) + wi(t - r3)). 
4 



where ei = Xlm=i "mttcT' ^2 = 

Now, the parameters Cj are linearly independent since 



^4 m 

=1 Ctm^l , £3 



Zll=i «m&r and £4 







£2 




£3 









,1 T71 



.3 TT- 



"1 



«1 

-rl 



hi 62 

&3 



3 



n 



bi 



^3 / 



0:2 
03 

\ "4 / 



/"I \ 



M Miag(r/i,77i, 773,773) 



Q;2 
as 
\ 0:4 J 



(5.10) 



We now compute the real D3-equivariant A-versal unfolding. Prom [2] we know that the real 
unfolding is given by 

C{a) = jCo + aiRe (Li) + Q;2lm (Li) + a^Re (L3) + a4lm (L3) 

so the delay-differential equation is identical to (5.9) with 

ei = ttiRe (oq) + Q;2lni (aj) + Q;3Re (oq) -|- Q;4lm (cq), 

62 — aiRe {a\) + a2lm. {a\) + a^Re (al) + a4lm [af) , 

63 — aiRc (1)2) + 02lin (bl) + a^Re (bl) + ailm (bl) , 

64 = ctiRe (63) + a2lra (63) + Q;3Re (6|) + a4lm (6|). 
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Moreover, the new parameters e^, i — 1, 2, 3, 4 are also linearly independent since 







( Re (aj) 


Im (aj) 


Re(ag) 
Re (a?) 


Im (og) 
Im (af ) 


\ 


/ «! \ 






Rc {a\) 


Im {a\) 




0:2 






Rc {b\) 


Im \b\) 


Rc {hi) 


Im 










\ {hi) 


Im {hi) 


Re 


Im (61) 


/ 


\ "4 / 



where the matrix is nonsingular since the four columns of the matrix span the four dimensional 
real subspace of the four dimensional complex space spanned by the columns of the matrix 
in (5.10). 

5.2 Double imaginary eigenvalues. 

We now consider the case where the eigenvalues of A are double. The space P has basis 

m = {Ue),UO),UO),UO).UO),UO),h{e),UO)) 

= {ve'^^\ ve'^^^, Ue-^^i^ ve''^^^ , ve'^^\ ve'^^^ , ve''^^^ , ve''^^^), 

where 




and UJ — e^™/^. We now compute the basis ^ of the adjoint P*. The action of D3 on is 
generated by matrices 

G{k) = diag (J, J, J, J) ' ' ' ' ' ' ' (5.11) 

where J = [iq]- Ftom equation (2.8), we know that ^ commutes with the actions (5.1) 
and (5.11). Now, recall that \& is chosen such that {^ , $) = / and note that ^3 — (pi, 04 = 
07 = 05, 08 = 06, then along with the equi variance condition on ^ we obtain 



*(0) = 



/Ml 


UJ 


UJ) 




Ml 


UJ 


UJ) 




71(1 


UJ 


UJ) 




71(1 


UJ 


UJ) 




C(i 


UJ 


UJ) 




C(i 


UJ 


UJ) 




C(i 


UJ 


UJ) 




V C(i 


UJ 


UJ) 


J 



where and ( are some nonzero complex numbers. 
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Basis of the unfolding space. Using (5.3), the dimension of the unfolding space W for 
matrix B is 16 since here r = 4 is the number of distinct eigenvalues, kj = 2 is the number of 
Jordan blocks for each eigenvalue and rij^e = 1 is the size of each Jordan block. 

We now define the 16 matrices which form the basis of S{yV). We define 3x8 matrices 
Rj;i,'&,m for each eigenvalue A^, where Ai = coii, A2 = — <^ii, A3 = uj2i, A4 = —u!2i- Let 
Vi — {vl,vf, vfY, i = l,...,8be3xl complex vectors. The linear mappings 11^ are 




(5.12) 



We choose ^,'(9 G {1,2} and m = 1 since each Jordan block is one-dimensional. Then, 
Rj;^,'&,m has vector V2{j-i)+^ in column 2{j — 1) + and all the other columns are zero where 




(5.13) 



and ^3 = vi, V4: = V2, vj = v^, vg = v^. 

To simplify the notation, we order the 16 matrices Rj;^,.fi,m in the following way: 

|^l+8(fe-l)^ ■ ■ ■ , = {Rl-l,k,l, R2-l,k,l, R3;l,k,l, R'i;l,k,l, -Rl;2,fe,l, R2;2,k,l, Rz;2,k,l, -R4;2,fe,l} 

for k — 1,2. One can verify using a computer algebra package that 

rank($(ro), $(ri), $(r2), $(t3)) = 8, 

therefore we need four distinct delays to,ti,T2, ts to solve 

3 

= ^A™$(ri). (5.14) 

1=0 

Projection of the basis. We project the matrices i?™ on the subspace of Ds-equivariant 
3x8 matrices, using the actions on (5.1) and (5.11). 

Proposition 5.1 Let v = {u^ , u"^ , u^Y G and consider the 3x8 matrix Np having v as its 
p*'* column and all the other columns are zero. From (5.5), we obtain a projection operator 
from which we compute 

Moreover, let ?7i(z^) — + uv^ + lov^ , r\2{y) — ^ lov^ + Ujv^ and recall that v — (1, a;, cJ)*. 
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1. If p — 1,5, then Np has only nonzero columns p and p + 1 given respectively by rjiv and 

2. if p — 4, 8, then Np has only nonzero columns p — 1 and p given respectively by vjiv and 

3. if p — 3, 7, then Np has only nonzero columns p and p + 1 given respectively by rj2V and 

V2V, 

4- and if p = 2, 6, then Np has only nonzero columns p — 1 and p given respectively by r)2V 
and r]2V. 

Proof The proof is a straightforward calculation. I 

Prom the definition of the matrices i?"*, for m — 1, . . . , 8 we have -& — 1 and for m = 
9, . . . , 16 we have § — 2. Table 1 shows the correspondence of matrices EJ^ with matrices Np of 
Proposition 5.1 as well as the information given by (5.12) and (5.13) needed to compute R . 
Prom Proposition 5.1 and Table 1, we obtain easily that R =0form = 5,...,12 and 

R^ = R^^ = [ji-^v, n-^v, 0, 0, 0, 0, 0, 0] , R^ ^ R^^ = [0, 0, ]J-^v, JJ-^v, 0, 0, 0, 0] , 

= = [0, 0, 0, 0, c-^^, C-'^, 0,0], '^^r'' = [o,o,o,oAO,r\ Vv] . 

Ds-mini-versal unfolding Since S{yV) is an unfolding space for B, we again set B^^ = 

^(0)i? for j = 1,...,16. The space W is spanned by sixteen 8x8 matrices. Por z = 1, . . . 8, 
fli has a 1 at element {i,i) and zeroes elsewhere. Por i = 9, 10, 11, 12, fli has a 1 at element 
{i — 8, 13 — i) and zeroes elsewhere while for i — 13, 14, 15, 16, J7j has a 1 at element {i — 8, 21 — i) 
and zeroes elsewhere. Now, obviously, B^^ = for j = 5, . . . , 12 and a computation shows that 

[B, yi] + (1 + (1 + i^-')\uj\^)n, + (1 + 2\uj\^)n2 + /U-V(l + 2|u;n(l^ii + O12) 
[B, 1/2] + (1 + 2|a;nn3 + (1 + 2\uj\^)n^ + {iJ4l-\l + \u;\^) + JI-'\uj\^)ng 
+yu7Z-^(l + 2|a;|2)fiio 

[B, ys] + (1 + (1 + C'')l^n^5 + (1 + 2\uj\^)n, + C-iC(l + 2|a;p)(l]i5 + Oig) 
[B, 1/4] + (1 + 2\u\^)Q-r + (1 + 2|a;|2)f]8 + (CC"'(1 + kl') + C~ Vl')^i3 
+CC~'(l + 2|cu|2)J^i4. 



r>D3 _ rDs _ 

-"1 — -'^13 ~ 

rBs _ dEDs _ 

-°2 — -"14 — 

hO-i _ rDs _ 

-^4 — -°16 ~ 
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r^nT'T'p*5"nn"nrlp'npp "wHfVi /V 


Pro 1 ppf i on 


7J - 




'/2 y^j) 


Rl /Vi anH /?9 ^ /Vo 








n 


r2 ^ aTq and /?^° •(-^ A/^^ 






n 




R3 ^. AT „„ J Dll ^ ]Vf 








n 


i?^ ^ A^7 and R^'^ ^ 


n4 


^7 







^ Ni and R^^ ^ N2 


Hi 


^^2 





/^-^ 


R^ ^ 7V3 and R^^ ^ N4 


n2 


V4 







R^ ^ ATg and R^^ ^ Nq 


na 


v& 







R^ ^ N7 and R^^ ^ Ns 


n4 




r 






Table 1: The first column shows the correspondence between i?™" and Np of Proposition 5.1 
while the remaining columns detail the information needed to compute matrices R of the first 
column where rij{-) is defined in Proposition 5.1. 

for some matrices yi,y2,y3,y4 G Matgxs- The four first rows of © are 

v^ = {l + {l + /i-i)|cc;|2,l + 2|a;|2, 0,0, 0,0, 0,0, 0,0,^-^(1 + 2|a;n,;U-i7l(l + 2|a;|2), 0,0, 0,0) 
V2 = (0,0,l + 2|cj|2,l + 2|a;|2,0,0,0,0,7I-i(/i(l + \Lof) + la;^), ^77-1(1 + 2|a;|2), 0, 0, 0, 0, 0, 0) 
V3 = (0, 0, 0, 0, 1 + (1 + 1 + 2|cj|2, 0, 0, 0, 0, 0, 0, 0, 0, ^^((1 + C^C(1 + 2|c^P)) 

V4 = (0,0,0,0,0,0,l + 2|cup,l + 2|a;|2,0,0,0,0,C"'(C(l + Ic^n + |cup), Cr'(l + 2|cj|'), 0, 0), 

rows 5 to 12 are zero and rows 13 to 16 are identical to rows 1 to 4. Thus, {fi, ^2, '^3, '^4} 
determine a maximal set of linearly independent vectors in the rowspace of ©. Hence, 

Mat^^g = TbTP' e C • {Bf^ Bf, Bf , B^^} 

and we need only consider matrices R' for j = 1, 2, 3, 4 in order to compute the Da-equi variant 
A-mini-versal unfolding. 

Computation of the matrices AJ*. We take Da-equivariant matrices as in (5.6) and 
we compute 
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Thus, equations 'W = X)5=o ^('O) 

m = 1, . . . , 4 can be solved by finding the solution to 

the first row only. 



A^^{Tj) for m = 


( 









= M 







V ) 





a\ + b\{uj 
al + blip 



UJ 



CO 



(5.15) 



where M is the generically nonsingular matrix (5.8). 

Remcirk 5.2 Note that, as in the previous example, from the form of equations (5.15) we can 
choose — a], = a^, — 6 - and b^ — b, for j — 0, 1, 2, 3. 



3-equivariant A-unfolding. For m = 


1,2,3,4, define 






3 


3 


( < w 


bf 


\ 




= E 


bf 


bf 




3=0 


j=Q 


[bf bf 


af 


J 



z{rj), 



and then the complex Da-equivariant A-versal unfolding is defined by: 



j=0 



( 



m=l 



b- 



bf 



bf 



bf \ \ 



z{rj). 



As in the previous example, setting tq = 0, ti = r*, T2 



T* and b'f^ 



= for m = 1,2,3,4 preserves the structure of the delay-differential equation. Then 
as a delay-differential equation, the complex unfolding is given exactly by (5.9) with again 



m=l 



a. 



e2 = Em=i «ma^, ea = Em=i OimK, and 64 = 



-I '-*m^3 



Note that the unfolding parameters ei, £2, 63, €4 are linearly independent by Remark 5.2 and 
from the same calculation as in the previous example. Moreover, from Remark 5.2, the real 
Da-equivariant A-versal unfolding is also given by (5.9). 

Appendix 

A Proof of Lemma 3.3 

We will prove the result in the case where (3.2) holds with p equal to the codimension of Tb^^ 
in Mat^x^. The proof follows along the same lines as the non-equivariant case (which can be 
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found, for example in [1]). We first define the F-centrafizer oiu E Mat^^^ as the hnear subspace 
of all matrices in Mat^^^ which commute with u. Consider now the F-centralizer Z£ of the 
matrix B. In the space GL(c, C) fl Mat^^^, let V he a smooth surface transversal to Jcxc + 
at /cxc, and such that the dimension of V is equal to the codimension of Z^ in Mat^^c- 
then define a mapping 

T-.TxCP — ^Mat^^^, J^{p,a) =pB{a)p-\ 

The key result we need is 

Lemma A.l T is a local diffeomorphism in a neighborhood of{Icxc, «o)- 
Proof First, we note that the mapping 

/ : GL{c, C) n MatLe ^ c Mat^e , 

/(C) = CBC-^ 

is such that the derivative of / evaluated at I cxc is 

Dfilcxc) : Tj^xM<xc = Mat^,, TbTF C T^Mat^e = Mat^,, , 

Df{Icxc){u)^[u,B], 

which is onto TbE'". It follows from the above that the dimension of Z^^ is equal to the 
codimension of TsE'" in Mat^^c) which is equal to p. Also, the dimension of V is equal to the 
dimension of ■ 

Now, the mapping T is such that 

DpJ^{Icxc, (^o){u, a) = Df{Icxc){u) = [u, B] , D^Tilcxc, Oio){u, a) = D^B{aQ)a . 

By construction of we have that 

^cXc 

maps Tj^^^V isomorphically onto TbS^. Also, by the hypothesis of Lemma 3.3, DaB{ao) 
maps TaoC^ = isomorphically onto a space which is a direct sum complement of TbE'" in 
Mat^x^. It thus follows that D{F{Icxc, chq) is an isomorphism between vector spaces of dimension 
dim(Mat^^p), so is a local diffeomorphism by the inverse function theorem. I 

Now, let A{(3) be a g-paramctcr F- unfolding of B (with A{Po) — B). Define Hi and 112 as 
the projections of P x onto V and respectively. 
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For all (3 sufficiently close to /3o in C^, define 



c(/3) = ni^-i 



{Am- 



It follows that for all (3 sufficiently close to I3q 



A{I3) = m) = C{(3) B{m) 



-1 



which proves the lemma. 
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